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Abstract — A Multi-hop Control Network consists of a plant 
where the communication between sensors, actuators and 
computational units is supported by a (wireless) multi-hop 
communication network, and data flow is performed using 
scheduling and routing of sensing and actuation data. Given 
a SISO LTI plant, we will address the problem of co-designing 
a digital controller and the network parameters (scheduling 
and routing) in order to guarantee stability and maximize a 
performance metric on the transient response to a step input, 
with constraints on the control effort, on the output overshoot 
and on the bandwidth of the communication channel. We 
show that the above optimization problem is a polynomial 
optimization problem, which is generally NP-hard. We provide 
sufficient conditions on the network topology, scheduling and 
routing such that it is computationally feasible, namely such 
that it reduces to a convex optimization problem. 

I. Introduction 

Wireless networked control systems are spatially dis- 
tributed control systems where the communication between 
sensors, actuators, and computational units is supported by 
a shared wireless communication network. Control with 
wireless technologies typically involves multiple communi- 
cation hops for conveying information from sensors to the 
controller and from the controller to actuators. The main 
motivation for studying such systems is the emerging use 
of wireless technologies in control systems (see e.g. [13], 
[12] and references therein). The use of wireless Multi-hop 
Control Networks (MCNs) in industrial automation results 
in flexible architectures and generally reduces installation, 
debugging, diagnostic and maintenance costs with respect 
to wired networks. Although MCNs offer many advantages, 
their use for control is a challenge when one has to take 
into account the joint dynamics of the plant and of the com- 
munication protocol. Wide deployment of wireless industrial 
automation requires substantial progress in wireless trans- 
mission, networking and control, in order to provide formal 
models and verification/design methodologies for MCNs. In 
particular, co-design of the control algorithm and of the 
network configuration for a MCN requires addressing issues 
such as scheduling and routing using real communication 
protocols. 

Recently, a huge effort has been made in scientific 
research on Networked Control Systems (NCSs), see 
[23], [7], [4], [15], [14], [11], [24], [19] and references 
therein for a general overview. In the current research on 
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the interaction between control networks and communication 
protocols, most efforts in the literature focus on scheduling 
message and sampling time assignment for sensors/actuators 
and controllers interconnected by wired common-bus net- 
works, e.g. [5], [8], [2], [22], and [18]. 

In general, the literature on NCSs addresses non-idealities 
(such as quantization errors, packets dropouts, variable sam- 
pling and delay and communication constraints) as aggre- 
gated network performance variables, losing irreversibly the 
dynamics introduced by scheduling and routing communi- 
cation protocols. What is needed for modeling and ana- 
lyzing control protocols on multi hop control networks is 
an integrated framework for analysing/co-designing network 
topology, scheduling, routing and control. In [3], a simu- 
lative environment of computer nodes and communication 
networks interacting with the continuous-time dynamics of 
the real world is presented. To the best of our knowledge, 
the only formal model of a MCN has been presented in [21], 
where the modeling and stability verification problem has 
been addressed for a MIMO LTI plant embedded in a MCN, 
with the controller already designed. A mathematical frame- 
work has been proposed, that allows modeling the MAC 
layer (communication scheduling) and the Network layer 
(routing) of the recently developed wireless industrial control 
protocols, such as WirelessHART (www . hartcomm2 . org) 
and ISA-100 (www.isa.org). The mathematical frame- 
work defined in [21] is compositional, namely it is possible 
to exploit compositional operators of automata to design 
scalable scheduling and routing for multiple control loops 
closed on the same multi-hop communication network. In 
[20], we extended the MCN definition to model redundancy 
in data communication (i.e. sending sensing and actuation 
data through multiple paths), with the aim of rendering the 
system fault tolerant with respect to link failures and robust 
with respect to packet losses. In [17] a MCN is defined as an 
autonomous system where the wireless network itself acts as 
a fully decentralized controller, and differs from the model 
we use in [20] and in this paper, where the wireless network 
transfers sensing and actuation data between a plant and a 
centralized controller. Moreover, in our model we explicitly 
take into account the effect of the scheduling ordering of 
the node transmissions in the sensing and actuation data 
relay, which provides more information when modeling the 
transient response of the system. 

In this paper we consider a MCN that consists of a 
controllable and observable continuous-time LTI SISO plant 
interconnected to a digital controller via two (wireless) multi- 
hop communication networks, as illustrated in Figure Q] 




by an arbitrary acyclic radio connectivity graph. In partic- 
ular, it allows modeling wireless multi-hop communication 
networks that implement protocols such as WirelessHART 
and ISA- 100. 



Fig. 1. Multi-hop Control Network scheme. 

An example of such networked control scheme is the mine 
application investigated in [6], where an industrial plant is 
connected to a controller via a wireless multi-hop commu- 
nication network. Since we assume that the data transfer is 
performed using multiple paths and merging redundant data 
by linear combination, the communication networks will be 
characterized by LTI SISO dynamics that depend on the 
network configuration (topology, scheduling and routing). In 
[20] we showed that, no matter how complex is the topology 
of the network, it is always possible to design a controller 
and a network configuration (scheduling and routing) that 
guarantees asymptotic stability of the networked closed loop 
system. In this paper we address the problem of co-designing 
a digital controller and the network parameters (scheduling 
and routing) in order to guarantee stability and maximize a 
performance metric on the transient response to a step input. 
More precisely, we will minimize the £2 norm of the error 
signal while requiring the closed loop system to follow a step 
reference with zero steady state error in finite time (deadbeat 
control), with constraints on the maximum overshoot in the 
input and output signal of the plant and on the maximum 
data rate on the communication links. 

We first assume that topology and scheduling are given 
and show that the above optimization is a polynomial op- 
timization problem, which is generally NP-hard, and we 
provide sufficient conditions on the network parameters such 
that it reduces to a convex optimization problem. Then, 
we consider the scheduling as a design variable, and we 
exploit the above optimization problem to define a total 
ordering induced by the £2 norm of the error signal on the 
finite set of scheduling functions. This allows to compute 
the optimal scheduling with a combinatorial computational 
complexity with respect to the number of communication 
links of the network. In future extensions of this paper we 
aim at providing algorithms to reduce such computational 
complexity. 

The paper is organized as follows. In Section|II]we provide 
the model of a MCN, while in Section [Til] we define the 
transmission data rate of each node. In Section HVl we define 
the co-design optimization problem and provide the main 
results of the paper. 

II. Modeling of MCN 

To allow systematic methods for designing the commu- 
nication protocol configuration in a MCN, a mathemati- 
cal model of the effect of scheduling and routing on the 
control system is needed. Definition Q] allows modeling 
time-triggered communication protocols that specify TDMA, 
FDMA and/or CDMA access to a shared communication 
resource, for a set of communication nodes interconnected 
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Fig. 2. Time-slotted structure of frames. 

In these standards the access to the shared communication 
channel is specified as follows: time is divided into slots 
of fixed duration A, and groups of IT time slots are called 
frames of duration T = IIA (see Figure [2). For each frame, 
a communication scheduling allows each node to transmit 
data only in a specified time slot. The scheduling is periodic 
with period IT, i.e. it is repeated in all frames. 

Definition 1: A SISO Multi-hop Control Network is a 
tuple Af = (V, Gn, Vn, Go, Vo, A) where: 

• V(s) is a transfer function that models the dynamics of 
a continuous-time SISO LTI system. 

> G-jz = (Vji, E-ji,Wti) is the controllability radio con- 
nectivity acyclic graph, where the vertices correspond 
to the nodes of the network, and an edge from v to 
v 1 means that v' can receive messages transmitted by 
v through the wireless communication link (v, v'). We 
denote v c the special node of Vjz that corresponds to 
the controller, and v u £ Vr the special node that 
corresponds to the actuator of the input u of V. The 
weight function W-r : Er — > R associates to each link 
a real constant. The role of W-r will be clear in the 
following definition of rjR. 

• rjn : N — > 2 Btc is the controllability scheduling func- 
tion, that associates to each time slot of each frame 
a set of edges of the controllability radio connec- 
tivity graph. Since in this paper we only consider 
a periodic scheduling that is repeated in all frames, 
we define the controllability scheduling function by 
rjiz : {1, . . . , IT} — > 2 Elz . The integer constant IT is the 
period of the controllability scheduling. The semantics 
of rjn is that (v,v') £ rj(h) if and only if at time 
slot h of each frame the data content of the node v 
is transmitted to the node v'. We assume that each link 
can be scheduled only one time for each frame. This 
does not lead to loss of generality, since it is always 
possible to obtain an equivalent model that satisfies this 
constraint by appropriately splitting the nodes of the 
graph, as already illustrated in the memory slot graph 
definition of [21]. 

« Go = {Vo, Eq, Wo) is the observability radio connec- 
tivity acyclic graph, and is defined similarly to Gr_. We 
denote with v c the special node of Vo that corresponds 
to the controller, and v y £ Vo the special node that 
corresponds to the sensor of the output y of V . 

• rjo : {1, . . . , IT} — > 2 E ° is the observability scheduling 
function, and is defined similarly to rj-ji. We remark that 
IT is the same period as the controllability scheduling 
period. 



• A is the time slot duration. As a consequence, T = IIA 

is the frame duration. 
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Fig. 3. MCN closed loop control scheme. 

For any given radio connectivity graph that models the 
communication range of each node, designing a scheduling 
function induces a communication scheduling (namely the 
time slot when each node is allowed to transmit) and a 
routing (namely the set of paths that convey data from 
the input to the output of the connectivity graph) of the 
communication protocol. Since the scheduling function is 
periodic the induced communication scheduling is periodic, 
and the induced routing is static. The above definition allows 
modeling redundancy in data communication, i.e. sending 
the same control data through multiple paths in the same 
frame. This kind of redundancy is called multi-path routing 
(or flooding, in the communication scientific community), 
and aims at rendering the MCN robust with respect to node 
failures and packet losses. To define the semantic of the data- 
flow in the networks we assume that each communication 
node computes a linear combination of the data received 
from the incoming links according to a weight function (i.e. 
W-r or Wo), and forwards these data to the next links. Note 
that this forwarding protocol does not require the use of a 
reordering protocol (such as TCP) and of a timestamp in the 
packets to indicate how recent is the data. This reduces the 
complexity of the protocol and the bit length of each packet 
(and thus the data traffic in the network). 

As derived in [20] the dynamics of a MCN N can be 
modeled as in Figure [3] where each block is a discrete-time 
SISO LTI system with sampling time equal to the frame 
duration, characterized by the transfer functions G-n(z), 
P(z) and Go(z). Throughout the paper, given a transfer 
function F(z), we will denote by Np and Dp respectively 
its numerator and denominator. We derive the discrete-time 
transfer function representation of the plant input/output 
dynamics P(z) by discretizing V(s) with sampling time T: 

, . _J_ P > ( \ = bn-iz 11 - 1 + b n _ 2 z n - 2 + . . . + bp 
[Z) M(z) [Z) M{z)(z r + a r ^ 1 z r - 1 + ... + a )' 

where M(z) is the polynomial whose roots are the stable 
poles of P(z). Since P(z) is derived by discretizing a 
continuous-time transfer function, the degree r + deg(M(z)) 
of the denominator is equal to n. As illustrated in [20], the 
transfer function of the controllability graph G-jz(z) can be 
expressed as follows: 

E m(d)z D ^- d 

G n {z)= d ^ zDii , 7*(d)= W M> 

(1) 

where xn(d) is the set of paths of Gu characterized by delay 
d, Dr C N is the set of the |D^| different delays introduced 
by the paths of G-n, Dn = max{D7?.} is the maximum 



introduced delay, W-ji(p) is the product of the weights of 
all links that generate path p of G-r,. We define the vector 
In = [ln{d),d € T> n ] G M |Dtc| . The transfer function of 
the observability graph Go{z) is derived similarly. 

According to the scheduling function, different paths can 
be characterized by different delays. In [20] we illustrated 
that, in order to avoid issues related to zero-pole cancelations 
between the controllability and observability networks and 
the plant (which leads to loss of controllability and ob- 
servability), it is convenient to design a scheduling function 
such that all paths are characterized by the same delay, i.e. 
\D-n\ = |De>| = 1. Moreover, as illustrated in SectionlTVl this 
situation is also convenient since it reduces the computational 
complexity of the optimization problem addressed in this 
paper. However, the following example shows that designing 
such a scheduling can increase the frame duration and 
therefore the total time delay introduced by the network, and 
it is not necessarily the best solution in terms of system 
performance, as illustrated in Example [4] When multiple 
control loops are closed on the same MCN it is even more 

difficult to design such a scheduling, as illustrated in [6]. 
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Fig. 4. Controllability graphs of Example^and[3](a), and Example [5] (b). 

Example 1: Let us consider the controllability graph il- 
lustrated in Figure Ufa) and assume that, because of inter- 
ference in the shared wireless channel, simultaneous trans- 
mission is only allowed for the set of nodes {vi, v$} and 
{v2,V3,V4}. Let us suppose A = 10 ms, and consider the 
following controllability scheduling functions rfj^ and 77^: 

VTzi 1 ) — {('"li^a)) (v i,t>a)j 0»ii "•»)}, 

1)^(2) = {(v 2 ,v 5 ), (v 2 , V 7 ), (1)3,115), (v 3 , v 6 ), (v 3 , 117), (v 4 ,v 6 ). (u 4 , v 7 )}, 
»?TC.( 3 ) = {(«B,«t). (v e - v 7 )}. 

'V-R.C 1 ) — {("li "2)1 ("1) °3)j "4), (V5, V 7 ), (l> 6 , V 7 )}, 

Vjii 2 ) — {(«2, us), ("2, v 7 ), (v 3 , K5), (u 3 , v 6 ),(v 3 ,v 7 ), (v4.,v e ), (,vi, v 7 )}. 

Using scheduling 77^, all paths are characterized by the 
same delay, namely G^(z) = -. Since IT a = 3, the total 
delay introduced by the network is given by IPA = 30 ms. 
Using scheduling 77^, paths are characterized by different 
delays, namely G h n (z) = °- 6z + 0A . Since n h = 2, the 
total delay introduced by the network is given by II^A = 
20 ms. It is clear that there is a non-trivial tradeoff in the 
design of the scheduling function to obtain the minimum 
total network delay while obtaining simple transfer functions 
for the controllability graph, which motivates our modeling 
framework. 

III. Data rate model of communication links 

Consider a control scheme where actuation (resp. sensing) 
data generated by the controller (reps, by the sensor node that 
measures the output of the plant) is quantized and transmitted 
through the controllability (resp. observability) network. We 



assume that such quantization process is characterized by 
a quantization width S u (resp S y ) over a range of values 
bounded by [-U max , U max ] (resp. [-Y max ,Y max ]). This is 
due to an implementation constraint, since the number of bits 
that we can use to represent the data is finite. 

Remark 1: In this paper, we only consider the effect of 
quantization as limited bandwidth constraints. It is out of the 
scope of this paper to characterize the effect of quantization 
on the dynamical properties of the closed loop system. For 
such results, the reader is referred to the existing literature on 
quantized control of linear systems, e.g. [16] and references 
therein. 

We assume that quantized data can be transmitted in the 
communication networks according to one of the following 
computational models implemented by the nodes. In the first 
model, each node computes the sum of all data received 
by the incoming links, and then transmits these data to the 
outgoing links multiplied by the corresponding link weight. 
In the second model, each node computes the weighted sum 
of all data received by the incoming links, and then transmits 
these data to the outgoing links. The advantage of the first 
model is that it allows each node to transmit with different 
data rates on the outgoing links, thus exploiting all available 
bandwidth on all links, while the second model restricts 
nodes to transmit at the same data rate on all outgoing links. 
On the other hand, as we will see in Section IIVI the second 
model induces weaker conditions for solving efficiently the 
optimization problem of the closed loop system performance. 
It is worth to remark that the above computational models, 
even if they affect in a different way the data rate modeling, 
generate the same controllability and observability networks 
dynamics as in (Q~|). 

Computational model 1: We first provide an example to 
show that the number of bits needed by a generic link (v,v') 
to transmit data depends on the weight function Wr. 

Example 2: Let us consider the controllability graph as in 
Figure|4jb) and suppose that 6 U = 1 and U max = 2. Suppose 
that W(vi,V3) = 1: the link (vi,vs) needs just 2 bits for 
representing the data. Suppose that W(vi,V2) — 10~ 3 : to 
avoid loss of precision when storing data in node v 2 we need 
to transmit on («i,i>2) more bits with respect to (vi,vs). 

Given a controllability network Gr = (Vr, Er., Wr), 
we propose a model for relating the transmission data rate 
of each link with the weight function Wr, and with the 
quantization parameters 8 U , U max , used to quantize the data 
at the input of the controllability network. For each link 
e = (v,v') G Er we define a variable a e as the sum of 
weights of all paths starting form v c and terminating with 
the link e: 



data divided by the transmission time A: 



d£D T . 



<X( v ,v')(d) 



53 Wr( P ), (2) 



p — v u , . . . , v, v 



The data rate R e associated to data transmission on link e 
is given by the number of bits needed to represent received 



R P = A" 



log 2 



2U n 



min \a e 



(3) 



where [•] is the ceiling function. Therefore, the data rate 
depends on the weight function Wr via equation (f2]). Note 
that also the variables / y 1z as defined in (Q]i are related to the 
variables a e as defined in 



53 a e (d), 



(4) 



e£inc(v u ) 

where inc(v) = {(v,v) G Er : v G Vr} is the set of 
incoming links in v. It is easy to check that the minimum 
achievable data rate is obtained for a e = 1, Ve € Er, which 
corresponds to the weight function defined by W(v, v') = 

\inc(v)\ ' w ') e Th' s cn °i ce i s not necessarily the one 
that optimizes the performance of the closed loop system, as 
illustrated in Section ITVl 

Computational model 2: Using the same reasoning as 
above, we define for each v G Vr \ v u a variable a v as 
follows: 



a v — 53 a v(d), a v (d) = 



d£D T . 



E 

p — v u , . . . , i 



Wr( P ), (5) 



with a Vc = 1 since the controller node directly handles the 
quantized data at the input of the controllability network. 
Furthermore, since the actuator node is directly connected 
to the actuator and therefore does not transmit data in the 
network, we do not define a„ u . The transmission data rate 
of node v and the variables jR(d),d G are given as 
follows: 



R v = A" 



log 2 



2U„ 



7n(d) 



v£pre(v u ) 



Su 

W(v,v u )a v (d), 



mm 



(6) 



(7) 



where pre(v) — {v : (v, v) G Er} is the set of prede- 
cessor nodes of v. It is easy to check that the minimum 
achievable data rate is obtained for a v = 1, Vi> G Vr \ 
v u , which corresponds to the weight function defined by 

E W(v', v) = 1, Vl> G Vr \ v u . 

v' Ginc(v) 

Depending on the use of the first or second computational 
model, we define the vector olr respectively by [a e ,e G 
Er] G M |£k| and [a v ,v eVr\ v u ] G Rl^l" 1 . The same 
reasoning can be used to characterize the data rates of the 
observability network. 

IV. Co-design optimization problem 

We consider in this section the problem of co-designing 
a digital controller C(z) and the weight functions Wr and 
Wo, in order to minimize the quadratic C 2 norm of the 
error signal e(k) = y(k) — r(k) while requiring the closed 
loop system to follow a step reference with zero steady state 
error in finite time (deadbeat control), with constraints on 



the maximum overshoot in the input and output signal of the 
plant and on the maximum data rate on the communication 
links. 

A. Controller and weight functions co-design for a given 
scheduling function 

The following example motivates our investigation by il- 
lustrating that, for a given scheduling function, minimization 
of the quadratic £2 norm of the error signal, of the overshoot 
and of the data rate on the links are conflicting specifications. 

Example 3: Let us consider a plant and a controllability 
graph as in Figure Ufa), described respectively by the transfer 
functions P(z) = ^ and G n (z) = 2&{affi&gl , i n 
this example we assume that communication nodes use the 
computational model 2 and, for simplicity, we will not 
consider the observability network. Consider a controller of 
the following form: 



C[z) = 



z 2 (d\z + do) 



(z - l)(z 2 + c\z + c ) 

To design a deadbeat controller we impose Np + Dp = 
z 4 , where F(z) — C(z)Gn(z)P(z), thus obtaining the 
following overshoot of the plant output: 

207 K (1)7r(2)+217k(1) 



O y = max 



3 7tc (1) 2 -7k(2) 2 +27k(1)7k(2) 
37k(1) 2 +207k(2) 2 +207tc(1)7k(2) 
37-r(1) 2 -7-r(2) 2 +27-r(1)7-k(2) 



Assume that S u = 0.1, U max = 500, A = 10ms and 
7tc(1) = I* then the overshoot is given by: 

20 7k (2) 2 + 20 7k (2) + 3 



O v 



-7 TC (2) 2 + 2 7K (2) + 3 



When a cancelation occurs with an unstable pole of the plant 



(i.e. 



7tc(2) _ 
7tc(1) 



3), then the overshoot tends to infinity. We 



discuss now two co-design cases: 
Case 1: Suppose to design W-r, so that 7 k(1) = 450 
and 7k (2) = 0.003. Designing the controller variables 
Co, ci, do, di to obtain deadbeat control, we obtain: 

II e(Jfe) ||£ 3 = 6.245, O y = 7, R V2 = 3.1kHz. 

Case 2: Suppose to design Wr., so that 7 tj.(1) = 1 and 
772.(2) — 1. Designing the controller variables Co, c\, do, d\ 
to obtain deadbeat control, we obtain: 



e(k) \\c 2 = 13.5508, 



Oy = 502.5, 



R V2 = lAkHz. 



The co-design of case 2, with respect to case 1, improves 
the transmission data rate of node V2 at the cost of strongly 
increasing the overshoot and the £2 norm. 

The above example shows that the co-design of controller 
and weight functions is non-trivial, and therefore motivates 
addressing the following optimization problems, for the 
computational models described in Section Hill 

Problem 1: Given a MCN TV, let V, (V K ,E K ), rfr, 
(Vo,Eq), rjo and A be given. Design a controller C(z) 
and controllability and observability weight functions Wji 
and Wo to minimize the quadratic £2 norm of the error 



signal of the closed loop system, subject to the following 
constraints: 

1) the closed loop system is asymptotically stable, with 
zero steady state error in finite time with respect to a 
step reference (deadbeat control); 

2) the overshoots of the signals u{k) and y(k) are norm 
bounded; 

3) the transmission data rate of each link (resp. of each 
node) is upper bounded. 

Since we consider a step reference input, we design a 
digital controller C(z) in the following form: 



C(z) = M{z)z D ^z Do C'{z), 

n t, x d s z s + d s - 1 z s ~ 1 

O (z) 



d 



(z-l)(z™ + c m _ 1 z™- 1 + ... + c ) 



(8) 



where M(z), z D71 and z Dtc are used to delete the stable 
poles of the plant and of the controllability and observability 
networks. In order to obtain the smallest response time of the 
closed loop system, we design C(z) with same degree in the 
numerator and denominator polynomials, therefore m + 1 = 
s + deg(M (z)) + Dn + Do ■ Note that the closed loop system 
is forced to follow a step reference with zero steady state 
error in finite time if and only if 



D C D P/ +N c ,N G ^NpN Gc 



(9) 



where / = m + r + 1 — Do is the response time. Moreover, 
we consider the weight functions W-jz and Wo as additional 
variables of our optimization problem. In other words, we 
address co-design of the controller C(z) and of the dynamics 
G-jz(z),Go(z) due to the data flow through the networks. 
The variables of our optimization problem are defined by 
the following real vectors: 

c = [c m _i,c m _ 2 , • • • , co]; d = [d a ,d a -i, ■ ■ ■ ,d ]; 
wr = [Wrc(e),Ve £ E n ]; w = [W©(e),Ve G E ]- 

Computational model 1: Problem Q] can be formalized as 
follows: 

min \\e(k)\\l 2 

c,d, w-re .wo 

subject to: 

1) D c Dp, + N c ,N G ^N P N Go = z l ; 

2) U < U , Oy < Oy, 

3) R e < R e , Ve s En U Eo- 

The following propositions characterize the algebraic relation 
between the objective function and constraints in Problem Q] 
and the optimization variables c, d, w^j, and wq. 

Proposition 1: The quadratic £2 norm of the error signal 
of the closed loop system is a polynomial function of d and 



Proof: The Z transform of y(k) in the closed loop 
system can be written as 

C(z)G n (z)P(z) 



Y(z) = 



l + C(z)G n (z)P(z)G (z) 
z D °N c ,N G7l N P 



D C D P , +N c >N G ^N P N Gi 
qi- v z l ~ u + . . . + qiz + q Q 



R(z) 
-R(z) 



R(z) 



(10) 



where v = min{D7^.} + 1 > 2 is the difference between 
the number of poles and the number of zeros of the closed 
loop transfer function, and R{z) — A-^j, A £ K is a step 
reference input. It follows that 

r, \ (qi-v , q\ Qa\ Az 



Y( 



qi 



l 



and 



(ID 



y(k) = Aqi- U 5^{k -u) + ... + Aq S-x(k - I), 
where <5_i(fc) is the unitary step. Therefore 
e(fc) = Aq^ v 8-i(k -u) + ... + Aq S-x(k -I)- AS-xQt). 
Thus, 



|e(*0| 



-A, 



-A,Aqi- v - A, Aqi_ v - 



-Aqi- v -x ~ A, . . . , ^ {Mh - A) 



h=0 



(12) 



As shown in ( [Tol l, the parameters q^ are sums of products 
of the coefficients of Nc, N G7l , Np, Dq as follows: 



qh = ^2 h ■ dj ■ jn(n), 

i+3 + (D^-K)=h 



(13) 



Proof: By (fTTT i the overshoot O y can be written as 
follows: 



Oy = 



-A,...,-A,Aq 



l-v 



A,Aq^ v + Aq, 



l-u-l 



l-v 

■,J2(Aq h 



A 



h=X 



A. 



(14) 



Since parameters qh are as in (13[ , then the constraint 

1 constraints | — A\ < 

E(Aq h -A) 



Oy < Oy can be rewritten as I — v 



O y ,\Aqi- v - A\ < O y ,.. 



h=l 



< O v , which 



are polynomial with respect to d and w^. Using the same 
reasoning, it is easy to prove that the same holds for the 
constraint on O u . ■ 
Proposition 4: Constraint 3 of Problem Q] is a polynomial 
constraint with respect to xv-r and we>. 

Proof: The constraint R e < R e , e £ E-r can be written 

as: 



Re = A" 



log 2 



2U„ 



mm Cte 



(2U max ^ \(X e | 
S u min \a s \ 



< 2 



\a e \ < 



S u 2 



< Re 



Re A 



— mm \a s 



K e min |a g | < \a e 



where -fn( K ), K £ T>k are polynomial functions of wjj as 
in (03. This concludes the proof. 

■ 

Remark 2: Note that, as shown in (fl~2l . the quadratic £2 
norm of the error signal increases as v = min{D7j} + 1 
increases. This implies that the minimum delay min{D7?.} 
is the one that mainly affects the responsiveness of the 
system to variations of the reference signal, and is thus more 
significant than the other delays introduced by the network. 

Proposition 2: Constraint 1 of Problem Q] is a polynomial 
constraint with respect to c, d, w-ji and wg. 

Proof: Equation (01 can be rewritten as follows: 

z l + ( P i- 1 -l)z<- 1 + ...+ Po = z l , 

where the parameters ph are given by 

Ph= ^2 Mi7R(«)7o(0+ a i + °j' 

i+J + (DK-K) + (Do-t.)=h i+j=h 

and where 77^. (k), k £ D-r. and 70 (t), t G Do are polyno- 
mial functions of w-ji and wq as in (Q~|). This completes the 
proof. ■ 
Proposition 3: Constraint 2 of Problem Q] is a polynomial 
constraint with respect to d and w^. 



K e \a s \ < 0, Ve e E n , 
(15) 

where n e > 0, and by (0 the variables a e ,e £ En are 
polynomial functions of w^. The proof for e £ Eq is 
similar. ■ 
Propositions Q]|4] show that Problem Q] is a polynomial 
optimization problem with respect to the variables c, d, 
w-r. and wo, which is well known to be NP-hard (see [1] 
and references therein). Therefore, unless P — NP, there 
exist no polynomial-time algorithm to solve any instance of 
Problem Q] In the following we provide conditions such that 
Problem Q] is a convex optimization problem. 

Lemma 5: Using the computational model 1, the follow- 
ing statements hold: 

For any given vector a-ji, there always exists a 
unique vector w-jj that satisfies equation (fJJ; 
Each component of 7^ is a linear combination of 
the components of a-ji if and only if 

Vei, e 2 <E inc(v), v £ Vr \ v u , 3\d E D n : 

Vp e {J xn(d), e x £ p, e 2 i p, (16) 

namely if and only if links belonging to paths 
characterized by different delays never merge in the 
same node. 
Proof: (i): By equation (0 it follows that: 



(1) 
(11) 



W n (v,v') 



(17) 



Since the graph is acyclic, for any given att ViV i) we can 
assign the variable Wn(v, v') so that equation (O is satisfied, 



without affecting the value of the variables auj >v \, v G 
inc(v), which can thus be considered as constants in equation 
(fTTT i. By applying this reasoning for all (v,v') G En the 
result follows. 

Sufficiency of (i): If (fT&t holds then each link (v,v u ) G 
E-R will belong to paths characterized by the same delay 
d v , namely at v<Vu \ = a/ ViVu )(d v ), This implies that j(d) = 

v:d v —d 

Necessity of ( ii): If condition ( fT6b is false, then there exist 
at least two links (vi,v') and (v 2 , v') that respectively belong 
to paths pi G Xn{d\) and p 2 £ Xn(d 2 ), di,d 2 G D TC . Then 
there exists at least one link (u, u„) € £7^ that belongs both 
to pi and p 2 - Let W-r,(V,u) = ^ Wk(p), then: 

a (v,v u ) = (a Vl {di) + a V2 (d 2 ))Wn(v' ,v)Wn(v 1 v u ), 
m(di) = a Vl (d 1 )Wn(v',v)W' R (v,v u ). 

Since a Vl = a Vl (di) and a V2 — a Vl (d 2 ), then jn(di) = 
. This concludes the proof. ■ 
The above proposition similarly holds for olq, -y and wp. 
A naive way to efficiently solve Problem Q] is by removing 
the dynamics on the networks: 

Proposition 6: Let \\)n\ = 1 and |Do| = L then Problem 
Q] is a convex optimization problem with respect to the 
variables c and d. 

Proof: If \D-jz\ = 1 and |Do| = 1, by Propositions 
[T]|4]it follows that the absence of j n and j decouples the 
optimization variables c, d and the variables wn and wo of 
the data rate constraints. By statement (i) of Proposition ( TToT ). 
it follows that the data rates can be arbitrarily assigned by 
designing appropriately wr and wq. Finally, the objective 
function and the constraints 1 and 2 of Problem Q] become 
convex with respect to c and d. ■ 
However, as we illustrated in Example [TJ the case \T>n\ = 1 
and |De>| = 1 is not easily achievable. For this reason, we 
consider the following relaxation of Problem [T] 

Problem 2: Let (d*,a^j) be the Pareto optimal value of 
the following optimization problem: 

min ||e(fc)||| 2 

subject to: 

2) O u < O u , O y < O y ; 

3a) R e < R e , Ve G E n . 
Given d = d* and an — a n , design c and olq such that 
the following hold: 

1) D C D P , + N C 'N G ^N P N Go = z l ; 

3b) R e < i? e ,Ve G E Q . 
If Problem [2] is feasible, then Problem Q] is feasible as well, 
and the Pareto optimal values of Problem [TJ and of the first 
part of Problem [2] are equal. This is due to the fact that the 
variables c and ~f do not appear in the objective function. 
Using the first computational model, we cannot state a priori 
conditions that guarantee feasibility of the second part of 
Problem[2] we will show that this is possible using the second 



computational model. If instead Problem [2] is not feasible, 
then Problem [T] might be feasible, and the Pareto optimal 
value of the first part of Problem [2] is a lower bound for the 
Pareto optimal value of Problem [TJ (if they exist). This is due 
to the fact that, splitting the Problem [T] we are neglecting 
some feasible solutions that require the joint design of all 
optimization variables. 

Proposition 7: Let w-^ > 0, wp > 0, let ( TTol l hold both 
for Gn and Go, and let |d| = 1 (namely the numerator 
of the controller is a polynomial of zero degree). Using the 
computational model 1, the first part of Problem|2]is a convex 
optimization problem with respect to a.n, and the second 
part is a system of linear equations with linear constraints, 
with respect to c and ocq. 

Proof: Since d| — 1 and each component of ~f n is a 
sum of components of ot-ji by Lemma [5] then the quadratic 
C 2 norm of the error signal and constraint 2 are convex with 
respect to ol-ji. Since w-ji > 0, wq > 0, then constraint 3a is 
linear constraint with respect to a-ji. Finally, since d = d* 
and otji — a^, constraint 1 is a system of linear equations 
in the variables c and olq, and constraint 3 is linear with 
respect to olq. ■ 
The above proposition is based on the assumption that 
wr > 0, wo > 0. If on one hand this makes Problem Q] 
convex, on the other hand it implies that 7^ > 0, j > 0, 
which reduces the set of designable dynamics for G-ji(z) 
and Go(z). Therefore, we might need to use more complex 
controllers and observability networks in order to make 
Problem Q] feasible. It is also assumed that links belonging 
to paths characterized by different delays never merge in 
the same node, which restricts the set of admitted network 
topologies and scheduling functions. We will show that, if 
the communication nodes use the computational model 2, we 
do not need to raise these two assumptions. 

Computational model 2: Problems [TJ and [2] can be 
formalized as above, by replacing the date rate constraint 
with R v < R v ,Vv G V n U Vq \ {v u , v c }. 

Lemma 8: Using the computational model 2, and for any 
given -f n and a^, it is always possible to design wr such 
that equations (f5J) and (0 are satisfied. 

Proof: It is easy to derive from (0 that, for all v G 
Vn \ {v c ,v u }, a v = W(v,v)ay. Since the graph 

vGpre(v) 

is acyclic, and reasoning as in Lemma [5] it follows that 
given ct-jz it is possible to assign a v , Vv G Vn \ {v c ,v u } 
by designing W(v',v) for just one v' G pre(v). All the 
remaining W(v, v), v G pre(v)\v' can be considered as free 
variables, that we will use to design -f n . Let pre(v u ) = F: 
note that the weights W(v, v u ), v G F have not been used to 
assign the vector an, therefore they can also be considered 
as free variables to design 7^. 

We claim that it is always possible to assign the vector 7 7^. 
using the free variables defined above. If |F| > \T>n\, men it 
is clearly possible to assign -f n by arbitrarily choosing any 
\T>n\ weights W(v,v u ), v G F. 

Assume instead that |F| < \Dn\- Let be the set of 
delays d G Dn such that there exists at least one node 



v E F with a v = a v (d), i.e. one node that belongs to 
paths characterized by a unique delay. By definition of 
it follows that |F| > |D^| + 1, since all paths associated to 
delays in do not merge before reaching node v u . Let 
Fi C F be the set of nodes belonging to paths characterized 
by a unique delay: it follows that we can assign jTz(d) for all 
d G by designing the weight Wn{v, v u ), for an arbitrary 
v € Fi such that a v = a v (d). 

Let D 2 n = D n \ D^, and let F 2 = F \ Fi. Any node 
v £ F2 belongs to paths characterized by a finite number 
M < |Dtj| of different delays. However, as discussed above, 
this implies that there exist at least /.i free variables associated 
to W(v,v u )a v , i.e. — 1 for all merges of different delays 
and one more for the variable Wr(v, v u ). By repeating this 
reasoning for all v S F2, it follows that we have enough 
free variables to assign "fn{d) for all d £ D 2 ^. ■ 
By the above lemma it follows that constraints 3 a and 3b are 
decoupled from Problem [2] and is therefore always possible 
to design wr and wo to assign the optimal data rates both 
in the controllability and observability graphs as illustrated 
in Section Hill 

Proposition 9: Let |d| = 1. Using the computational 
model 2, the first part of Problem [2] is a convex optimization 
problem with respect to 7 K , and the second part is a system 
of linear equations with respect to c and 7 . 

Proof: If |d| = 1, by Lemma [8] it directly follows that 
the quadratic £2 norm of the error signal and constraint 2 are 
convex with respect to -y n , and that constraint 1 is a system 
of linear equations in the variables c and -f . Moreover, by 
Lemma[8]it follows that constraints 3a and 3b are decoupled 
from the objective function and all constraints of Problem [2] 
and therefore they disappear. ■ 
The following proposition provides a sufficient condition on 
topology, scheduling and routing of the observability graph, 
such that Problem [T] and Problem [2] are equivalent. 

Proposition 10: Let |d| = 1 and \Dq\ > r + 1. Using 
the computational model 2, Problem [T] and Problem [2] are 
equivalent. 



Proof: The condition De> > 



1 guarantees the 



existence of a solution on the variables c and 7 for the 
system of linear equations defined by constraint 1. Since 
by Lemma [8] the constraints 3, 3a and 3b disappear, and 
since the variables c and -f only appear in constraint 1, 
then constraint 1 can be decoupled from Problem Q~| and can 
also be decoupled from the second part of Problem [2] As 
a consequence, Problem Q] and the first part of Problem [2] 
become equivalent. ■ 

B. Scheduling design 

Problem 3: Given a MCN N, let V, (Vr., E n ), (V ,E ) 
and A be given. Design scheduling functions 777^ and rjo, a 
controller C(z), and controllability and observability weight 
functions Wr, and Wo, to minimize the quadratic £2 norm 
of the error signal of the closed loop system, subject to the 
constraints of Problem Q] 

The set T of admissible scheduling functions is finite, and 
its cardinality is exponential with respect to the number of 



communication links and to the scheduling period. Therefore, 
Problem [3] can be solved by computing the Pareto optimal 
value £\{rfR, 770), ~i(r\R, rjo) S T, using the results obtained 
in the above section. This induces a total ordering in the set 
r, and we can define the set of optimal scheduling functions 

r * = {(Vn,V*o) e F: (Vr,V*o) = argmin £* 2 (m,Vo)}- 

(>jTC,>jo)er 

Example 4: Let us consider a plant, a controllability graph 
G-ji, and parameters S u , A and U max as in Example [3] We 
also consider here an observability graph Go equal to G-r, 
and the following observability scheduling function 

»)o(l) = {(vi, v 2 ), (vi, l> 3 ), (i>i, v&), (v 2 ,v s ), (v 2 ,v 7 ), 

(V3,v 5 ), (v 3 ,v 6 ), (1)3,1)7), (V4, v e ), (v 4 ,v 7 ), (v 5 ,v 7 ), (1)6,-07)}, 

that is associated to the following transfer function: 

7o(2)z + 7o(3) 
Go(z) = r; . 

We consider the following three scheduling functions: 

'JtcW — {(^i)^), V3), (1)1,1)4), (v5,vr), (ve, vr)}, 

ri^(2) = {(112,1)5), {i>2, V7), (v 3 ,v 5 ), (113, ve), (v 3 ,v 7 ), (v 4 ,v 6 ). (1)4,1)7)}, 

*?"£.(!) = {(fit V2), («1, D3), (till "4), (l)2, D7), (l) 3 , D7), (l)4, U7), (D5,D7), 

(ve, vr)}, 

77 TC (2) = {(1)2, U5), ("3, V5), (1)3, v e ), (D4, v b ), }, 

VtiW — {{vi, v 2 ), (i>i, D3 ), (" 1, 14), (1)2, 1)5), (f3,D5), ("3, v e ), (1)4, v e ), 

(1)4, 1)7), (1)6,1)7)}, 
1tz( 2 ) = {(v2, vr), (1)3,117), ("5, 1)7)}, 

that are associated to the following transfer functions: 



G K °(z) 
G nb (z) 
G n .{z) 



7K (l)z + 7R (2) 
z 2 

7 k(1)z 2 + 7tc(2)z + 7k (3) 



Let O u = O y = 10, and 1.5 kHz < R e < 3 kHz, V e £ 
Er U Eo- We chose a controller with |d| = 1 and |c| = 
4 for rj^ and 77^, and with |d| = 1 and |c| = 5 for 
rfj^. We assumed that the communication nodes use the 
computational model 1, and solved Problem [T] for each 
scheduling using CVX, a package for specifying and solving 
convex programs [9], [10]. 
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Fig. 5. |J e(k) ||£ norm vs. maximum allowed data rate, (a) solid line: 
rj^\ (b) solid line: 17^; (b) dashed line: rfj^. 



Figure [5] shows the £2 norm of the error as a function 
of the data rate bound R e , for the three scheduling defined 
above. Scheduling rjL is characterized by a single delay, and 
has the worst performance, probably because it introduces 
in the optimization problem only one degree of freedom. 
However, rf^ performs worse than 77^ despite the fact that 
it introduces more degrees of freedom, probably because the 
additional free variables are redundant to reach the minimum 
£2 norm and only increase the network delays leading to 
a worse performance. In conclusion, the scheduling that 
guarantees the best £2 norm while respecting the constraints 
is 77^. Figure [5] also shows that it is useless to increase the 
transmission data rate on the links to values greater than 
1700kHz, since the performance does not improve after that 
point. 

V. Conclusions 

Given a MCN, we addressed the problem of co-designing 
a digital controller and the network parameters (topology, 
scheduling and routing) to guarantee stability and maximize a 
performance metric on the transient response to a step input, 
with constraints on the bandwidth. We showed that the above 
problem is a polynomial optimization, which is generally 
NP-hard, and we provided sufficient conditions such that it 
reduces to a convex optimization. 
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